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ABSTRACT 

This  paper  considers  the  time  dependent  Stefan  problem  with  convection  in 
the  fluid  phase  governed  by  the  Stokes  equation,  and  with  adherence  of  the 
fluid  on  the  lateral  boundaries.  The  existence  of  a  weak  solution  is  obtained 
via  the  introduction  of  a  temperature  dependent  penalty  term  in  the  fluid  flow 
equation,  together  with  the  application  of  various  compactness  arguments. 
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SIGNIFICANCE  AND  EXPLANATION 

Consider  a  phenomenon  (such  as  melting  of  ice)  where  there  is  a  change  of 
phase,  say  liquid-solid.  In  the  liquid  phase  the  thermal  energy  is 
transported  both  by  diffusion  and  convection,  and  the  effects  of  convection 
are  reflected  in  the  movement  of  the  free-boundary  separating  the  two  phases. 

y1  In  ihis  paper  we  show;  that  such  a  problem  can  be  formulated 

A 

mathematically  and  that  it  admits  a  solution  in  a  weak  sense. 

11  "c" 

y — tie  faleo  investigate^some  local  regularity  properties  of  the  distribution 
of  temperature  and  the  field  of  velocities  in  the  liquid  phase. 


The  responsibiity  for  the  wording  and  views  expressed  in  this  descriptive 
summary  lies  with  MFC,  and  not  with  the  authors  of  this  report. 
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THE  BIDXMENSXONAL  STEFAN  PROBLEM  WITH  CONVECTION i  THE  TIME-DEPENDENT  CASE 


J.  R.  Cannon11*,  e.  DiBenedetto12*  and  G.  H.  Knightly13* 

1.  INTRODUCTION 

Th«  ala  of  thia  pa par  ia  to  axtand  to  the  nonatationary  caaa  some  raaulta  obtained  by 
tha  authora  in  [5,  6] ,  about  tha  Stefan  problem  with  convection. 

He  briefly  deacribe  tha  phyaica  of  tha  phenomenon,  referring  to  aection  2  for  a 
praciaa  mathematical  formulation. 

Suppoae  that  in  a  region  Q  of  RM,  N  >  2  a  liquid  undergoee  a  change  of  phaae  at  a 
fixed  temperature.  The  model  example  we  have  in  mind  ia  a  water-ice  aituation.  At  every 
time  t  the  liquid  and  aolid  phaae a  are  determined  by  the  Knowledge  of  the  diatribution  of 
temperature,  we  call  u1 1  *  tha  temperature  in  the  liquid  and  u12*  the  temperature  in 
tha  aolid.  in  general,  in  the  liquid  region  there  are  preaent  convective  motiona 
originated  by  body  forcea  I  depending  on  the  temperature  u1 1  * .  The  dynamic  atate  of  the 
liquid  ia  determined  by  the  knowledge  of  the  field  of  velocitiea  v  and  the  praaaure  p. 
The  diffuaion  of  heat  in  the  liquid  ia  affected  by  the  velocity  v,  and  in  turn  v  point 
by  point  ia  affected  by  the  buoyancy  forcea  ?(u* 1  * ) . 

He  will  deacribe  the  phenoawnon  of  diffuaion  in  the  liquid  phaae  by  the  evolution 
equation 

(1.1)  Vu"*)  -  div  k1Cu<1’)Vxu<1)  ♦  v-V^Cu11*)  -  0 

where  cf(*),  and  ( • )  repreaent  heat  capacity  and  conductivity  reepectively  and  are 
poeaibly  nonlinear  functiona  of  the  temperature  u11 *.  The  term  v*V  qj  (u1 1 ' )  gives  a 
deacription  of  how  the  velocity  v  affecta  the  temperature  u11 * . 


1 1 *The  University  of  Texaa  at  Auatin,  Texas,  78712 
12 'Indiana  Univeraity,  Bloomington,  Indiana,  47405 
13 'university  of  Masaachuaatta  at  Amherat,  Massachusetts,  01003 

Sponsored  by  the  United  States  Army  under  Contract  No.  DAAG29-80-C-0041 . 


Th«  convection  will  be  modeled  by  the  system  of  Stokes  equations 


(1.2) 


it  *  ■  ^  *  V  *(u<1>)  ' 


.<1>, 


where  v  is  the  kineaatic  viscosity,  p  the  pressure  and  f(u  )  the  buoyancy  forces. 

The  two  equations  (1.1)-(1.2)  represent  the  classical  Boussinnesq  coupling  of  thermal 
diffusion  and  convection  (15). 

He  assume  the  liquid  is  incompressible  (div  v  ■  0).  Moreover  since  it  is  viscous 
(v  >  0)  and  since  we  assume  that  the  solid  phase  is  at  rest,  it  is  reasonable  to  assume 
v  “  0  on  the  boundary  of  the  liquid  region. 

In  the  solid  phase  there  is  only  a  diffusion  process  described  by  an  equation  like 
(1.1)  without  the  tern  involving  the  velocity,  since  we  assume  zero  velocity  for  the  solid 
phase. 

He  assume  the  distribution  of  temperatures  u*1*,  u*2*  and  the  field  of  velocities 
are  known  at  soon  initial  time  t  »  0,  and  on  the  boundary  3Q  of  Q  we  prescribe  at 
every  time  t  the  heat  flux  g,  which  is  a  possibly  nonlinear  function  of  the 
temperature. 


At  the  unknown  boundary  r  separating  the  two  phases  we  inpose  the  relation 


(1) 


(2) 


0  and 


[k1(u(1))Vxu(1>-k2(u(2>)7xu(2>].Nx-I,Nt 


where  N  =  (Nx<Nt)  is  the  unit  normal  to  T  directed  toward  the  solid  phase.  Such  a 
relation  measures,  roughly  speaking,  the  amount  of  heat  used  in  the  melting  process  and 
L  represents  the  latent  heat  of  fusion. 

The  problem  consists  in  determining  at  every  time  t,  the  distribution  of 
temperatures,  the  field  of  velocities,  the  pressure  and  the  configuration  of  the  system. 

Our  purpose  is  to  show  that  such  a  problem  for  the  spatial  dimension  H  -  2  admits  at 
least  a  solution,  in  a  sense  to  be  made  precise  below. 

He  comment  here  on  the  restriction  N  “  2,  and  on  the  difficulties  of  extending  the 


results  of  [5,  6] ,  to  the  time-dependent  situation. 


Since  the  Stoke*  equation*  have  to  hold  only  in  tha  liquid  region,  on*  ha*  to  have 
some  topological  information  on  tha  set  occupied  by  the  fluid,  in  order  to  give  a 
meaningful  interpretation  to  the  field  of  velocities .  For  example  on*  should  know  that 
such  a  set  is  open  to  view  the  Stoke*  equation*  at  least  in  the  sen**  of  distribution*  over 
such  a  sat.  This  information  would  be  implied  by  the  continuity  of  the  temperature  which 
in  turn  depends  on  the  smoothness  of  v.  It  turns  out  that  only  for  N  ■  2  are  we  able  to 
show  that  the  degree  of  smoothness  of  v,  suffices  to  yield  the  continuity  of  u. 

This  delicate  interplay  between  the  regularity  of  u  and  the  regularity  of  v,  has 
also  prevented  us  from  using  the  Navier-Stokes  equations 

v  -  vdv  (v*V)v  ♦  -  f 

in  the  place  of  (1.2), 

The  limitations  are  quit*  clearly  of  technical  nature  and  it  is  our  hop*  to  remove 
them  in  the  future. 

Also  it  should  be  pointed  out  that  at  this  stage,  uniqueness  is  an  open  question. 
Sections  2  and  3  contain  the  classical  formulation  of  the  problem,  notation  and  some 
preliminary  material.  The  concept  of  weak  solution  is  introduced  in  Section  4,  whereas 
Section  S  is  devoted  to  a  listing  of  the  assumptions  and  a  statement  of  the  results.  In 
Section  6  we  prove  our  theorem  by  assuming  certain  facts  (Propositions  6.1,  6.2,),  which 
are  demonstrated  in  Sections  7,  8. 

It  is  a  pleasure  to  acknowledge  conversations  with  Prof.  B.  Benjamin  and  w.  Pritchard, 
on  the  physic*  of  the  problem. 

2.  FORMULATION  OF  THE  PROBLEM 

2 

Let  ft  be  a  bounded  domain  in  R  with  smooth  boundary  3ft.  For  all  t  «  [0,T], 

T  >  0  let  0<t)  2  ft  x  {t>,  3fl(t)  -  3fl  x  {t}  and  0  -  U  0(i).  Me  denote  with  30. 

'  0<T«t  T 

the  parabolic  boundary  of  0^  i.e. 

3ft_  =  ft«f)  U  s„  ,  S  «  U  30(t)  . 

o<t«r 


The  set  B  is  divided  into  Q.  and  a.  by  the  fra*  boundary  r  £  T  £  u  Tit), 

T  12  T  0<t<r 

whara  P(t)  is  a  hyparsurface  in  Q(t)  determined  implicitly  by  4(x,t)  “  0.  The 
function  •  d  c1(Qt),  •  <  0  on  B^ ,  ♦  >  0  in  and  |  V^*|  *  0  on  f.  Hara  ?x 

danotas  tha  gradient  operator  with  respect  to  the  space  variables  x  £  (x^x^J  only. 

Tha  sat  HO)  divides  tha  initial  region  0(0)  into  two  regions  Qf(0)  and  8^(0). 
We  sat  SA  ■  n  and  denote  with  tha  outward  normals  to  S^. 

Consider  tha  problem  of  determining  the  real  valued  functions  0,  u^  i  8^  ♦  R,  a 
*  2 

vector  valued  v  t  0^  R  ,  and  p  i  8^  ♦  R,  satisfying 

(2.1)  a1<u,i))  -  div  ki(u(i>)Vxu<1)  +  v*Vi(u(i>)  “  0  iB  fli'  1  “  1,2 


(2.2)  -k1(u(’*’)  ’M*11  -  g1(x.t,u(i)),  ( x,  t )  eSj,  i-  1,2 


uli)(x,0)  -  u^Nxl,  x«  0i(0)  ,  (-dSi^Nx)  <  0 


"o  VcO)  “  ° 


lk1(u(1))Vxu<1>  -  k2(u(2))Vu<2,)*Vxo-  h*t,  (x,t)  e  r  , 


v  -  V4v  +  -  £(u{1))  in  a, 


v  •  0  on  S. 


v(x,0)  -  VQ(X),  X  e  a^O),  div  VQ  »  0 


vQ(x)  £  o,  x  e  Bjto) 


div  v  ■  0 


v(x,t)  “  0  a. a.  in  Bj  . 


Here  the  latent  haat  L  la  a  known  poaitlva  constant,  g^x.t.O,  1  ■  1,2  are  known 

functions  of  thair  arguments  and  £(  •)  la  a  given  vactor  valuad  function  napping  R  into 

2  i 

R  .  also  ( • ) ,  ki(*)  arc  anooth  functions  dafinad  for  (-1)  a  <0  and  such  that  thara 

axista  constants  Yfl,  Y,  for  which 


(2.10)  0  <  Y0  <  a'(s),  k^s)  <  Y,,  i  -  1,2 

V  s  «  «,  (-1)ls  <  0  . 

In  ordar  to  fornulata  tha  problan  in  a  sinpla  fashion  wa  aaka  a  change  of  tha  unknowns 
as  follows.  First  note  that  by  thair  physical  nature  tha  haat  capacities  or  ( ■ )  are 
continuous,  nonotona  increasing  and  coercive  functions  of  thair  arguments  which  can  be 
defined  so  that  0^(0)  >0,  i  ■  1,2.  Therefore  tha  knowledge  of  a^(u^)  determines 
u<il  uniquely.  Than  wa  define 


-1(u(,)) 


«2(u(a)) 


in  0. 


in  0. 


K(u) 


/  k,(«,  (C))«,  (5)dt  ,  u  >  0 


I  k .(a"1  «))<*"  <{>d«  ,  u  <  0 

0  * 


u0(x) 


g(x,t,u) 


<u’(x>) 

in  0,(0) 

(u£2,(x)) 

in  02<0) 

:,t,a''(u)) 

(x,t)  8  S 

:,t,a"V(u)) 

(x,t)  8  S 

-5- 


The  aquation*  for  u*1^,  i  -  1,2  can  be  formally  rewritten  as 

d  ♦ 

U  -  AK(u)  +  V*V^U  “0  in 

»  g(x,t,u),  (x,t)  8  ST 


(2.11) 


u(x,0)  -  uQ(x),  x  8  0(0) 
u(x,t)  »  o  (x,t)  e  r 

UVxk(u))+  -  tvxic(u)]_}*vx*  -  L*t,  (x,t)  e  r 

where  [V^KIu)]*  denotes  the  limit  from  0^,  as  (x,t)  approaches  r,  while  t V^Klu) ) 
denotes  the  limit  from  02* 

As  for  the  velocities,  setting  ?(u)  _  £(a^(u>),  we  can  rewrite  (2.5)  as 

(2.5)'  v  -  vAv  +  V^p  -  f(u)  in  0<  . 

In  order  to  formulate  our  notion  of  weak  solution  of  (2.11),  (2.5)',  (2.6)-(2.9),  we 
need  to  introduce  some  basic  notation. 


3.  NOTATION  AMD  FUNCTION  SPACES 

In  this  section  we  give  a  brief  description  of  the  function  spaces  employed  and  recall 

basic  facts,  known  from  the  literature,  to  be  used  as  we  proceed. 

For  q,r  >  1  let  L  (0  )  denote  the  Banach  space  of  those  measurable  functions 
q,  r  T 

napping  0T  ♦  K  with  norm  defined  by 


(3.1) 


lul 


where 


-6- 


,U,2eA<T)  -/  • 


When  q  -  r  »  2,  L  (0  )  coincides  with  the  Hilbert  space  L  (Q )  whose  inner  prod  :t 
2#  2  T  2  T 


(♦,•) 


,  „  generates  the  nora  1*1  „  «  1*1,  ,  „  .  11  <9  *  *  *• 

•  *>UT  i»2*uT 


set 


lul  n  ■  lul  „  • 

q.q.a,  q»aT 


Let  Wp'°(0T)  denote,  for  p  >  1,  the  Banach  space  with  nora 


(3.2) 


lulp  „  -  lulp  o  f  IV  ulp  -  , 

p-3t  x  P'V 

p  T 


where 


'Vp.Q,  */  IVIP<»* 


t  a 


and  | * |  here  denotes  the  euclidean  length  of  a  vector  in  R2. 

If  p  •  2,  is  the  Hilbert  space  with  inner  product 


(3.3) 


(u,w) 


Q  -  (u,w)  a  ♦  I  |j-)2 

2'°T  i-1  *i  *X1  2' 


2  T 


With  we  denote  the  Hilbert  space  with  inner  product 


(3.4) 


(u,w) 


<u,w) 


f 9u  9wi 

♦  ITT'  77 


. *;'H>  lw  ^2'°x 


3ti  du 

Here  denote  generalised  derivatives.  These  definitions  are  modified  in  the  usual 

way  if  q,r,p  are  infinity.  For  p  >  2  let  vT'0^),  be  the  subspace  of  ' 0 ( )  of 
those  functions  who**  trace  on  2Q(t)  vanishes  for  a.e.  t  9  [0,T1.  Also  let  p(AT) 
denote  the  Banach  space  with  nora 


(3.5) 


,U,V  (0  )  "  ,up  ,u,2.0(t)  *  ‘Vp.a  * 
2,p  t'  0<t<T  ' '  ’  P'  T 


-7- 


-8- 
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Corollary  3.2:  Let  Q  CO.  and  u  6  V.  (Q),  p  >  1,  Then  for  every  cylindrical  doaain 

T  2,p 

such  that  C  5 

,U,2P(Q-  *  C,U'v2<pttt) 
where  C  depends  upon  p  and  dist(Q,Q' ). 

Proof:  Q  has  the  fora  Q  =  a  »  where  G  is  a  region  in  Q  and 

0  <  tj  <  tj  <  T.  Let  K  be  a  compact  contained  in  G  and  set  Q'  <  K  *  (t^.t^). 
Construct  a  smooth  cutoff  function  9  CQ(G)  such  that  f  1  1  on  K.  Then 


u^j  9  V  (Q)  and 
2:P 


^  2P»Q  ^  ^  V,  <Q> 

2*P 


The  corollary  follows  from  the  particular  construction  of  *. 

* 

Next  we  describe  the  function  spaces  where  v  will  be  found. 

Denote  with  J(Q )  the  closure  in  the  nora  of  tL^ffl) ] 2  of  Vg{Sl),  the  space  of 

infinitely  differentiable,  divergence  free  vector  fields  $,  coapactly  supported  in  8, 

Let  also  J, (8 )  denote  the  closure  of  in  the  nora 

i  o 

Set  J(8t)  -  L2(0,TiJ(Q))»  W  -  L2<0,T|JJ(a>),  J10(8T)  •  C(0,T)J(8>)  ^  0,(8^, 

J1,liaT)  “  lW2,1(aT)l2  °  W  VV  “  n  VV*  «°te  that 


j“caT)  c  [V2(0T)12,  and  J,  0<aTJ  c  {v’'°<nT))2. 


4.  THE  WEAK  FORMULATION 

Let  9  Wj'  (8^,)  8UC**  that  v(x,T)  »  0.  Then  multiplying  the  first  of  (2.11) 

by  v  and  performing  (formal)  Integration  by  parts,  we  obtain 


//  {-8(u)Vt  ♦  V^K(u)*Vx*  +  v'V^u^Jdxdt 


-/  g(x,t,u)^do  +  /  S(u. )*(x,0Jdx  , 
ST  8(0)  0 


I  *, 


where  0  ( • )  is  the  eaximal  eo  no  tone  graph 


(4.2) 


0  ( s ) 


s  ,  s  >  0 

(— L,  0)  ,  s  ■  0 

s  -  L  ,  s  <  0 


The  formal  calculations  leading  to  (4.1)  are  routine  and  we  refer  to  [7,  12],  for 
details.  Here  we  only  remark  that  the  juep  of  0(‘)  at  zero  takes  into  account  the 
interface  relation  in  (2.11). 

Since  the  graph  0(*)  is  nulti valued,  0(u(x,t))  has  to  be  interpreted  <is  a  function 
w(x, t)  C  0(u(x,t)),  the  inclusion  being  intended  in  the  sense  of  the  graphs.  In  order  to 
siaplify  the  syabolisa  we  will  keep  the  synbol  0(u(x,t)),  bearing  in  Bind  the  way  it  has 
to  be  interpreted. 

Since  uQ  *  0  except  on  T,  0(uQ(x))  is  unambigously  a.e.  defined  in  Q(0). 

To  obtain  a  weak  formulation  of  (2.S)',  (2.6)-(2.9),  consider  a  smooth,  divergence 
free  vector  valued  function  $  which  is  coapactly  supported  in  Q^t)  for  all  t  6  [0,T] 
and  t(x,T)  S  0.  Take  the  “dot*  product  of  (2.5)'  by  $  and  integrate  by  parts  in  . 
Routine  calculations  [9,  10],  give 


<4-3)  //  {-v*J  +  vV “  f(u)*$}dxdt  -  /  v  (x)$(x,0)dx  . 

a,  x  0,(0)  ° 

Note  that  0^  is  the  set  where  u  >  0. 

Definition;  By  a  weak  solution  of  (2.11),  (2.5)',  (2.6)-(2.9),  we  aean  a  pair  (u,v)  such 
that 

(i)  u  8  v2(flT>  n  c(flT> 

(ii)  v  0  J1(0T)»  v  ■  0  a.e.  on  the  set  tu  <  0)  =  {(x,t)  6  0T|u(x,t)  <  0)  . 

(iii)  u  and  v  satisfy 


(4.4)  //  {-0(u)*  +  V  K(u)»V  *>  +  (v*V  u)^}dxdt  ■  -/  g(x,t,u)^do  ♦  /  0(u„ )l(x,0)dx 

«T  ST  0(0)  0 

for  all  *  0  w’*1  (0T)  such  that  i(x,T)  =  0,  and 


-10- 
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(4.5)  //  {-v-l  ♦  VV  V,V  {  -  £(u)$)dxdt  •  /  v  (x)|(x,0)dx 

(U>01  c  0,(0) 

for  all  <!/  e  J,  ,(Q_)  such  that  V  •$  -  o,  $(x,T)  •  0  and  tupp  $(»,t)  C  tu  >  0)(t). 
Remarks »  (1)  sinca  we  require  u  to  be  continuous,  the  set  [u  >  0]  is  open  in  the 

relative  topology  of  0^,  and  therefore  the  last  integral  identity  is  well  defined. 

(ii)  The  integrals  in  the  identity  of  the  teaperature  are  well  defined,  aodulo  some 
basic  assumptions  listed  below. 


5.  ASSUMPTIONS  AMP  STATEMENT  OF  RESULTS 

With  respect  to  the  data  g(x,t,u),  uQ(x),  vQ(x),  f(u)  we  assume  the  following 
[A,]  uQ(x)  *  0  •«*«  in  0(0),  Ug(x)  >0,  x  €  0,(0),  Ug(x)  <  0  on  02(0). 

Moreover  u0  is  essentially  bounded  and 

,uo'«,0(0)  *  * 0 

where  Kg  is  a  known  positive  constant. 

[a2]  The  function  g(x,t,£)  is  continuous  over  sT  *  ft,  and  satisfies  ths  growth 
condition 

|g(x,t,S)|  <  *„  +  K,ltl 

where  K,  is  a  known  positive  constant.  Moreover  t  ♦  g(x,t,C)  is  monotone  at  the  origin 
for  all  (x,t)  €  S  ,  i.e.  g(x,t,£)sign€  >  0. 

r  i 

[A3]  with  respect  to  K(')  we  assume  that  S  ♦  K(6)  is  Lipschits  continuous  in 
H\{ 0} ,  it  is  monotone  increasing,  and  satisfies 

0  <  XQ  <  K'(C)  <  X,,  a.e.  5  €  R\{0}  , 


where  XQ  -  and  X,-^*"1. 

IA41  vq(x)  €  J(Q)  and  vQ(x)  »  0  a.e.  in  OjlO). 

(a5]  f  is  Lipachitz  continuous  in  R  and 

|f(s,)  -  l(s2)|  <  Kjis,  ■  *2! 
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for  some  constant  K2  and  all  6  R,  i  -  1,2.  We  will  also  assume  f(0)  -  0.  This  is 
no  loss  of  generality  since  f(0)  -  0  can  be  always  be  realised  by  addition  of  a 
conservative  force  incorporated  in  the  pressure  tent. 

let  a  -  <<Va2)  €  (a  u  >  *  <■  u  f0*>*  be  a  multi-index  of  length  |a|  -  a,  +  «2* 
formally  with  D°F  we  denote  the  derivatives  of  r  of  the  form 


>1 


*1  °2 
3xl  ®X1 


P  . 


Also  for  o,n  €  (0,1)  let  H*^(°T)  denote  the  space  of  those  functions  which 

are  Holder  continuous  on  compact  subsets  of  0^  with  Hfflder  exponents  o  with  respect  to 

the  space  variables  and  n  with  respect  to  time. 

We  can  now  state  our  main  result 

Theorems  Under  assumptions  [A^J-IAj]  problem  (2.11),  (2.5)',  (2.6)-(2.9)  has  a  weak 
solution.  Moreover 


(i) 

u  a 

L„(Q_) 

■  T 

(ii) 

if 

5  ♦  K(5)  8  c” (R-)  then  i 

(iii) 

if 

e  ♦  K(C)  a  c“(B+J  and  if 

(a) 


d“u  a  i^oc(Bt  n  I u  >  OJ) 


Dau  e  H f'ff/2(0,  n  (u  >  01),  o  -  0(0)  e  (0,1)  . 
x  loc  T 


X  ♦  v(x, t)  e  c"(8  n  (u(* ,t)  >  0)),  a.e.  t  8  [0,T]  , 

(b) 

D°v  a  L_OC ( 0_  n  (u  >  OJ),  for  every  o  . 
x  •  T 

Remarks:  (i)  Assumptions  are  somewhat  stronger  than  what  is  needed  for  the 

proof  of  the  theorem  and  have  been  formulated  in  order  to  simplify  the  arguments.  In  fact 
the  monotonicity  of  C  ♦  g(x,t,€)  at  zero  can  be  relaxed,  and  on  u0(x)  one  only  needs  to 
assume  uQ(x)  8  i*2(Sl ( 0 ) ) .  We  will  indicate  later  how  this  can  be  done. 


(ii)  If  in  (2.11)  the  variational  boundary  data  ara  raplacad  by 

uL  -  h(x,t),  <*,t)  e  S 

ST 

than  we  have  a  Dirichlet  problem  whoaa  weak  foraulation  can  ba  darivad  in  an  analogous 
way.  However,  it  ia  nacaasary  to  taka  taat  functions  «  8  Xj*  Wj)*  W>*  proof  of  tha 
chaoraa  carrias  over  to  this  situation,  aodulo  tha  obvious  changes  due  to  tha  different 
nature  of  tha  boundary  data.  Me  cat  it  the  details  for  tha  Dirichlet  problaa. 

6.  PROOF  or  THg  TjMOggt 

Tha  plan  is  to  obtain  the  solution  u,v  as  a  limit  of  nets  (ut) .  (ve)  solutions  of 
certain  approximating  problems  solved  in  all  0^.  Since  v  must  act  ultimately  only  on 
the  sat  (u  >  OJ,  we  introduce  in  the  approximating  process  a  penalisation  acting,  roughly 
speaking,  on  (u  <  0]. 

Let  e  >  0  ba  fixed  and  consider  tha  problem  of  finding  a  pair  (ue,ve>, 
ue  a  W.  V  8  J*(0T)  satisfying 

II  {-0(u  )*  ♦  V  K(u  )•▼  P  ♦  v  »V  u  <P)dxdt  - 

-  -  /  g(x,t,ue)Wdo  I  0(uoMx,O)dx 

S_  0 

T 

*  8  «l'*(a_)  such  that  *(x,T)  5  0,  and 
2  T 

II  {-v  +■  W  v, »V  $  ♦  e'V  (u_>v,J  -  ?(u  >$}dxdt  -  }  v  !x)$(x,0)dx  , 

J.  C  t  X  *  X  *  *  *  *  Qw 

T 

^  €  Ji  il!V'  3  °*  H,re  e'V>  •  p*n*lty  t,tm  d*fino<1  ** 


(6.1) 

for  all 

(6.2) 

for  all 
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He(«) 


1  ,  -*■<•<  -X 

-e_1«  -  1  ,  -2e  <  •  <  -e 

0  f  -e  <  8  <  •  , 


♦  ♦  * 

«nd  vQ(x)  h«r«  is  th«  extension  of  vQ (x)  via  0  into  0^(0). 

The  crux  of  the  matter  lies  in  the  proof  of  the  following  facts. 

Proposition  6.1»  For  all  e  >  0  the  system  (6.1)-(6.2)  admits  a  solution  (ue,v£). 

moreover  there  exist  constants  co'C1'C2  dePend^n9  upon  the  data  in  assumptions  [Aj]-[Ag] 

but  independent  of  6  such  that 


(a)  *Vv  (0  )'  "V  - 

v2l  V  j,<a> 

1  T 


<  C. 


(b)  Iff  H  (u  )|£  |2dxdt  <  C. 
6  n  66  6  5 


(c)  'V-.aT  «  c2  ‘ 


For  every  compact  K  c  0^,  there  exist  constants  Cj,  C4  depending  upon  K  and  e , 
such  that 

ld)  ,ltV2.IC<  C3“C*‘> 
ta)  'V-.K  *  c4(K'e>  • 

Proposition  6. 2i  The  net  is  equicontinuous  over  0^ j  that  is  for  every  compact 

K  C  B^,  there  exists  a  nondecreasing,  continuous  function  «^(. )  •  R+  *  U+, 

«^(0)  »  0,  dependent  on  dist( tC.30^,)  but  not  upon  e,  such  that 

lUgU^tj)  -  ue(x2,t2)|  <  -  x2l  |tt  -  t2l/j> 

for  every  pair  (xi,ti>  «  K,  i  -  1,2. 

By  Proposition  6.1,  the  nets  {u  },  {v  }  are  weakly  compact  in  v*].'0 (ft  ) ,  and 

C  C  2  T 

J1 (0^)  respectively.  Hence  for  a  subnet  relabeled  with  e 
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uc  ♦  u  weakly  in  and 

v£  ♦  v  weakly  in  J^J^)  * 

By  proposition  6.1(c)  and  Proposition  6.2,  ths  nst  {ut}  is  aquibounded  and  equicontinuous 
on  ovary  compact  K  C5.  Thsrofors  a  subnet  can  bo  salsctod  and  rolabad  with  e  such 
that 

u  *  u  uniformly  on  compacts  K  C  Q_ 
c  T 

and  consequently 

u£  *  u  strongly  in  LjtO^) 

Because  of  the  equicontinuity  of  {u  }  the  uniform  limit  u  is  continuous  in  II  ,  and 

e  T 

therefore  the  set  (u  >  0]  is  open  in  the  relative  topology  of  Qf. 

The  proof  of  these  propositions  is  lengthy  and  technically  involved.  He  postpone  it 
to  the  next  sections  and  show  here  horn  to  conclude  the  proof  of  the  theorem,  using  these 
facts. 

[6.x] .  The  identity  of  the  temperature. 

Since  {u^}  is  aquibounded  in  0^,,  from  the  definition  (4.2)  of  the  graph  B,  it 
follows  that  the  net  (B(u£)}  is  also  aquibounded  in  B  .  Therefore  the  selection  of 
subnets  can  be  made  in  such  a  way  as  to  insure 

BlUj)  ♦  w  weakly  in  I.2(BT)  . 

By  aonotonicity  of  B(*)  we  have  w  c  B(u)  in  the  sense  of  the  graphs.  Xs  before  we  will 
write  B(u)  instead  of  w. 

By  the  trace  theorem  [14],  for  every  n  >  0  there  exists  a  constant  C(n)  such  that 

/  |ue  -  u|2do  <  nlV^(ue  -  #|lj  q  +  c(n)lue  -  ul2  . 

S  '  T  '  T 

T 

Since  the  norms  IV  (u  -  u)l  .  are  aquibounded  and  u„  ♦  u  strongly  in  L„(fl_),  the 
X  E  i,W  t  2  T 

T 

above  implies 

^  ♦  u  strongly  in  L j(S^) 
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and  consequently  in  view  of  the  continuity  of  g(x,t,u> 

g(x,t,ue)  ♦  g(x,t,u)  strongly  in  . 

By  assumption  [A3] ,  the  definition  of  K( •)  and  the  equicontinuity  of  {ue>  over  compact 
subsets  of  0^  it  follows  that  {K(ue>}  is  equicontinuous  on  compacts  K  C  and  hence 
the  selection  of  subnets  can  be  made  to  include 

K(Ug)  ♦  x(u)  uniformly  on  compacts  K  C  , 

Consequently 

K(ue>  ♦  K(u)  strongly  in  L2(£^)  * 

Now  we  also  have 

VK(ue)  ♦  s  weakly  in  L^J^)  • 

Let  e  C?(Q_).  Then 
0  T 

//  V  K(u  J*dxdt  -  -//  K(u  ) V  Abedt 


and  letting  e  ♦  0 

//  z*  -  -  //  K(u)  V  tf>dxdt  . 

°T  *T 

This  implies  s  =  V^Ktu).  For  the  nonlinear  term  involving  v£  in  (6.1)  we  have, 

//  v  *V  u  edxdt  “  -//  v  *V  #u  dxdt  ♦  -II  v*V  eudxdt  "  //  v*V  u*dxdt  , 

oT  ^  e  “  e  ^  x  ^  x 

for  all  yiew’'1(QT). 

Finally,  letting  c  ♦  0  in  (6.1)  we  obtain  the  identity  (4.4)  for  all  'f  e  ' 1  ( (^ ) 
such  that  l(x,T)  =  0. 


[6.B]  The  identity  of  the  velocity 

Let  K  be  a  coaipact  contained  in  the  open  set  [u  >  0).  Since  u(  ♦  u  uniformly 

on  K,  there  exists  eQ  >  0  so  small  that  ue(x,t>  >  0,  V(x,t)  «  )C  and  all  c  <  tQ. 

Let  $  6  *  <V,  $<x,t)  =  0  and  supp  ^  ( • ,  t )  C  K,  then  for  the  term  involving  the 

penalization  in  (6.2)  we  have 

e-1  //  Hc(u_.)v  l  dxdt  -  0  V  e  <  e.  . 

6  6  6  0 
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Letting  e  ♦  0  in  (6.2)  with  the  indicatad  choica  of  $,  and  Baking  uaa  of 
proposition  6.1,  wa  obtain  (4.5).  as  K  C  [u  >  01  waa  arbitrary,  wa  obtain  (4.5)  for  all 
auch  I  with  aupp  J(*,t)  C  (u  >  0](t). 

It  raaains  to  show  that  v  ■  0  a. a.  on  (u  <  0]. 

By  Propoaition  6.1 

jIJ  He(uc)|5e|2dxdt  <  c,  . 

°T 

Let  K  be  a  eoapact  containad  in  (u  <  0).  Proa  tha  unifora  convargance  of  u£  to 
u,  it  followa  that  thara  exists  so  saall  that  He(ue)  S  1  on  K  for  all  c  <  c^. 

Tharafore 

If  |v  |2dxdt  <  ec  . 

K 

Latting  e  ♦  0,  by  lower  seaicontinuity  of  the  nora  wa  obtain  v  -  0  a.e.  in  K,  and 
aince  (  C  |u  <  O]  ia  arbitrary,  v  «  0  a. a.  on  tu  <  0]. 

(6.0  Regularity 

Stateaent  (i)  will  be  proved  in  Proposition  6.1. 

(ii).  On  the  set  (u  <  0],  v  ■  0  a.e.  Therefore  the  tenperature  satisfies 

a 

u  -  &K(u)  -  0  , 

in  the  aense  of  distributions  over  [u  <  0).  Note  that  here  we  have  used  the  definition 
(4.2)  of  the  graph  B(»).  For  every  cylindrical  doaaln  Q  C  n  (u  <  0)  we  have 
u  9  Vj'°(Q)»  Therefore  if  *(•)  is  infinitely  differentiable  on  (-*,0),  the  stateaent 
is  a  consequence  of  classical  results  (ill. 

(lit).  On  the  set  (u  >  0],  u,  v  satisfy 

u  -  Ak(u)  ♦  v»7^u  ”  0 

a  ♦  ♦  _  ♦ 

v  -  viv  7^  -  f(u) 

in  tha  sansa  of  distributions  over  (u  >  0],  for  a  aaaaurable  function  p.  Moreover  for 

every  cylindrical  doaain  Q  C  (u  >  0)  we  have  u  6  V^'^tQ),  and  by  Proposition  (6.1), 

*  2 

u  fl  L  (Q)  and  v  8  IL  (fl)J  .  Hence  tha  stated  regularity  follows  froa  the  results  on  the 

•  © 
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local  smoothness  of  weak  aolutiona  for  tha  bidimensional  Bouaalnneaq  ayataai  established  in 


HI. 

7.  PROOF  or  PROPOSITION  6.1 

Denote  with  8  (•)  a  sequence  of  smooth  functions  in  R  such  that  8  (a)  ♦  8(a) 

Dfc  ■ 

over  compact  subsets  of  R\{0}  and  satisfying 
(i)  8  are  monotone  increasing 

(ii)  B^(s)  >1,  V  8  8  R 

(iii)  8  (s)  -  s  for  s  >  — ,  and  8  (a)  -  s  -  L  for  s  <  -  ^  , 
a  m  m  m 

Such  a  sequence  can  obviously  be  constructed. 

Let  also  (v  }  be  a  sequence  in  J.(Q)  such  that 

Ot  II  I 

*0,m  *  *0  in  J(2>  * 

Me  fix  e  >  0  and  for  each  m  8  ■  consider  the  following 
Auxiliary  Problem:  Find  u^  8  V^'0^).  *aeJi,i(flT>  s*tiBtying 


(7.1) 


//  {-ea<uB)*t  +  V(UB)*V  +  V'xV)dxdT 


-  /  g(x,t,u-W>do  +  /  8a(u0)'F(x,0)dx 


v  ^  8  wl'^Q-,)  such  that  «(x,T)  =  0 


2  T 


(7.2) 


//  {-  v  +  vV  v  >7  J  +  7  H  (u  )v  •$ 
„  mt  xmx  eemm 


f(u  )$}dxdT  -  /  v  *$(x,0)dx 

m  U  f  » 


V  $  e  J,  . (0  )  such  that  $(x,T)  =  0  . 
1  r  I  T 


Me  make  the  convention  of  denoting  with  C  a  generic  positive  constant  depending  upon 
quantities  that  will  be  specified  as  the  constant  appears. 
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Proposition  7.1:  For  each  m  6  M  the  auxiliary  problem  (7.1)-(7.2)  has  a  solution. 
Moreover  there  exist  constants  C  depending  upon  the  data  in  assumptions  [A^]-[Ag],  but 
independent  of  a  and  e ,  such  that 


(7.3) 


(7.4) 


lu  I..  Ivl  .  <  C 

*  ■  e  /  n  % 


-  vv 


W 


T  //  «  (u  )|v  |2dxdT  <  C  . 
el  e  m  m 


proof  of  Proposition  7. 1 i  He  employ  a  Galerkin  procedure.  In  L^tQ),  introduce  the 
orthononsal  basis  (z^(x)}  generated  by  the  problems 


-At^  ■  in  0 


\ 

in 


on  )Q 


In  J(Q)  we  introduce  the  orthonomal  basis  generated  by  the  Stokes  problems 


+  Vp  -  |i1Ci  in  0 


v*c1  -  0 


i  30 


where  P  is  a  scalar  function  representing  a  pressure.  Prom  (10,13)  it  follows  that 

form  a  complete  orthonormal  set  in  J(fl).  By  Kilbert-Schmidt  theorem  (13)  we  see 
that  any  smooth  divergence  free  vector  valued  function  n,  compactly  supported  in  0  has 
an  absolutely  and  uniform  convergent  representation 


n(x)  -  I  • 


i-1  2' 

Moreover  the  derivatives  of  n(x)  have  an  absolutely  and  uniformly  convergent 
representation  obtained  by  term  by  tern  differentiation. 
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i  ) 


He  represent  the  initial  data  u^fx),  vo,«^x*  as 


u0(x)  “  2  cixi(x) 


i-1 


"o,.(x)  ‘  j,  d&(x)  ‘ 


For  i  8  *  fixed  eet 


(7.5) 


v|(x,t)  -  l  di(t)Ci<x) 


and  denote  by  H  ,(x,t)  the  unique  eolution  of 
*#4 


(7.6) 


Jt  Vi  -  ♦  ?rv:1(Vi)  ■ 0 


(7.7) 


■VxK<8«1<%,i))*"  “  g<*'t'8«i1<,,n.i>>  °B  S,« 


(7.8) 


Vitx'°*  ”  vj*X*  ' 


in  the  eenae  of  projections  over  the  span  of  '*2'** * '*i^ '  where 


Vx)  “  l 


i-1 


and 


Wx”  ■  j,  • 


Namely  denoting  with  the  LjtB)  projection  onto  the  linear  span  of  {z, 

W  is  the  unique  element  in  N^IQ-)  satisfying  (7.8)  and 

*  2  T 


l  tit  VtV  +  W1(Vi,,,v*V  + 


(7.9) 


♦  v}*7  B"1  (N  ( )P,^}dxdt  -  -/  g(x,t,8"  (N  ,)?.«« 

i  x  a  m,i  i  '  m  «.*  * 

ST 
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1 


for  all  f  e  w2(a>,  and  all  t  e  [0,T]. 


For  tha  construction  of  W  .  (global  in  time) ,  we  refer  to  [2,3,4]. 

Me  will  employ  the  function  w  .  so  obtained  to  construct  the  function  v,  solution 


of 

(7.10) 


fe  *t  '  V^t  *  ■  f<0m  <Vi,)  in  *P 


(7.11) 

(7.12) 


vt(x,t)  “  o  (x,t)  e  sT 

v,(x,0)  -  l  d?c.(x)  , 

*  i-1  1  1 


in  the  sense  of  projections  over  the  span  of  .  Cj, . . . , Cj).  Namely  denoting  with 

the  J(Q)  projection  onto  the  linear  span  of  , Cj,...,  Cj),  v{  is  the  unique  element 

in  j  (8  )  satisfying  (7.12)  and 
1  #1  T 


(7.13) 


/  vt  nj  ♦  ,  vxntJ  ♦ 


for  all  $  e  (8)  and  for  all  t  e  [0,T]. 

Me  omit  the  proof  of  the  existence  and  uniqueness  of  satisfying  (7.10)  since  the 
construction  follows  by  straightforward  modification  of  standard  techniques  [10,9], 

Denote  by  Ls[0,T]  the  space  of  essentially  bounded  real  valued  functions  in 
[0,T],  and  with  H1 [0,T]  the  space  of  those  square  susmable  real  valued  functions 
t  ♦  d(t),  whose  weak  derivative  d'(t)  is  square  summabla  over  [0,T]. 

For  tel  set 

“  (Lw(0,T]  )*,  the  cartesian  product  of  t^lO.T) 
by  itself  t  times 

Yj  -  (H^O.T])*,  the  cartesian  product  of  H^tO.T] 
by  itself  t  times. 
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■  aai  am a  -  -■  - 


F 


In  Xg  and  Yg  introduce  the  norms 


id,  (•  ),d  (•),...  ,d,  (•  )l  »  ess  sup  £  fd.(t)]4 

12  *  Xt  0<t<T  i-1  1 


'VW’ . 'VllO.fl  \l  ,al'l  (0,T) 


The  procedure  described  above  defines  a  map  Fg  :  Xg  ♦  Yg  as  follows.  Given  an 

*  *  « 

l-tuple  ld1(t)«d  (t)>...ld((t)}  6  Xg  we  construct  g(x,t)  the  unique  solution  of 

♦  .see, 

(7.9)  with  v*  given  by  (7.5).  Then  we  associate  with  (d1 (t)»d2< t ), . . . ,dg(t )) ,  the 

Pourier  coefficients  {d  (t),d  (t),...,d^(t)}  of  the  solution  Vg(x,t)  of  (7.13). 

We  need  to  show  that  for  each  1  8  >,  F  possesses  a  fixed  point  in  Xg. 

This  can  be  done  by  using  the  Leray-Schauder  fixed  point  theorem  (10).  For  this  we 

have  to  prove  that 

(i)  Fg  maps  a  bounded  set  B  in  Xg  into  itself 

(ii)  Fg  i  B  ♦  B  is  compact 

(iii)  The  set  of  solutions  of  XPg(x)  ■  x,  X  6  (0,1)  is  bounded  independently  of 

X. 

Lemma  7. 1 :  There  exists  a  constant  C  which  depends  only  upon  uQ(x) ,  0^,  L,  the 
constants  in  assumptions  and  which  is  independent  of  e,  m,  1,  v*  such  that 

,Wm.t<*'t,,2  n  *  C'  a11  8  (0'T1  • 

Horeover  setting  0  1 (W  ,)  »  U  ,,  there  exists  a  constant  C  dependent  on  the  previous 

in  m,  *  in  t  ** 

♦ 

quantities  and  independent  of  e,  m,  1,  vj  such  that 


(7.14) 


♦J 

0 


IV  0  ,(• 

X  m 


,T)l2,ndT 


«  c>  t  e  (o,T) 


Proof  of  Lemsa  7.1;  The  proof  is  exactly  the  same  as  the  proof  of  Lemmas  1  and  2  of  [4]. 

♦ 

Here  we  only  show  why  the  constants  C  are  independent  of  v|.  The  lemma  is  proved  by 
setting  in  (7.9) 


V 


a,t 
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f 


and  performing  estimate*  in  the  identity  so  obtained.  For  the  term  involving  v*,  we  have 

l  *»£>.. t’Vm.i*'  ■  l  h{\  [n,t  C(5)d5)d*  • 

-  -/  [/  B,t  e*\t)de]div  ^*dx  -  o  . 

Q  0  * 


Lemma  7.2:  (a)  There  exists  a  constant  C  which  depends  upon  the  constants  in  Lemma  7.1, 

the  Lipschitc  constant  of  ?(•),  vQ(x)  and  which  is  independent  of  e,  m,  i,  v* , 

such  that 


(7.15) 


'V 


,VxVll*''t),2,Qdt  <C» 


t  e  [o,t] 


(7.16) 


e  //  H£<Um,t),^il2dxdT  *•  c 


(b)  Thera  exists  a  constant  C  which  depends  on  m,  e,  lvQ  ^  1^ 
independent  of  1  such  that 

i 

(7.17)  l|^  vtl2  <  c(m, e)  . 


but  which  la 


Proof  of  Le— n  7.2:  To  prove  (7.15)  and  (7.16)  choose  "  v^  in  (7.13).  To  prove 

A  f  ^ 

(7.17)  choose  v^.  The  lemma  now  follows  from  routine  calculations  and  Gronwall's 

inequality. 

Now  (7.15)  inqplies 

t 

l  (d  (t)] 2  <  C,  for  all  t  e  [0  ,T]  , 
i-1  1 

and  therefore  if  (d*(t) ,d*(t) , . . . ,d*(t) }  belongs  to  the  ball  of  radius  C^2  in  x^  we 
have 

Fj{d*,d*,...,d*>  e  the  ball  of  radius  C^2  in  . 

in  feet  by  (7.17)  ^  ,d2> . . .  ,dt)  -  Ft(d*,d*, . ..  ,d*}  e  Yt  and  Yt  is  compactly  embedded 
in  Xj.  Therefore  F^  is  compact.  The  continuity  of  Fj  is  demonstrated  by  a  standard 
difference  argument,  and  condition  (ill)  is  trivial  (see  (5)  for  similar  arguments). 


r 


Hence  we  conclude  that  F  t  actait*  a  fixed  point  in  Xt  and  actually  every  fixed 
point  liee  in  Y^.  Let  v(  be  i  fixed  point  of  ^  and  set 


Then  H  ,,  0  , ,  v, 

s,i  »,i  4 


“■.4 

satisfy  the  identities 


6.  * 


17.18) 

for  all 


<p  8  such  that 

2  T 


>v.* +  " 

*(X,T)  5  0 


-/ 

S 


g(x,t,0  , 
■»  — 
T 


)Pg*dc 


//  (|r  v.*H  $  +■  v7  v  i  7  I  J  ♦ 

LJ  l«t  4  s  xt  x  a 

mT 

(7.19) 

+  e"V(U  ,)v  n  J}dxdt  -  //  f(U  ,)*n  Jdxdt 

e  ■,<  x  S  Q  *1*  s 

T 

for  all  J.  ,  (Q  >  such  that  5(x,T)  3  0.  Here  s  is  a  positive  integer  <  4.  Now  we 

1.1  T 

let  4  ♦  *•  while  n  6  ■  remains  fixed. 

The  lieit  process  of  4  ♦  •  in  (7.18)  is  carried  out  exactly  as  in  [4].  He  remark 
that  a  crucial  fact  in  this  connection  Is  to  show  that 


(7.20) 


ug  strongly  in  L^O^) 

"e.4  *  w»  8  W'  W*“kly  in  W 


These  facts  were  shown  in  (4)  to  which  we  refer  for  details. 

The  passage  to  the  lieit  in  (7.19)  follows  the  argunents  of  [9.10],  by  making  use  of 
the  information  (7.20). 

Note  that  by  lower  senicontinuity  of  the  norm  the  estinates  (7.14),  (7.15),  (7.16), 
(7.17)  are  still  valid  in  the  lieit.  Proposition  7.1  is  proved. 

7- (a)  Regularity  of  u,,  and  v^ 

Next  we  give  an  equivalent  formulation  of  (7.1),  by  using  the  Steklov  averagings  of  a 
function  P  8  V2^T*'  defined  by 
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0  <  t  <  T  -  h 


€. 


P.  (x,t)  <• 
n 


t+h 

£  /  P(X,T)(Jt, 


0  t  >  T  -  h 

In  111]  it  ia  shown  that  (7.1)  is  equivalent  to 


(7.21)  //  ♦  VK(u-> VV  *  tVVxttsVldx<lT  ”  / 

°t  St 

for  all  0  <  t  <  T  -  h  , 
and 

8b(ub(x,0))  -  8b(uq(x))  , 

for  all  ♦e»!'t(l)J. 

2  T 

La— a  7.3»  Thera  exists  a  constant  C  independent  of  e,  such  that 

ass  sup  |u  |  <  C  m  »  1,2,...  . 

°T 

Proof «  Let  k  be  a  positive  real  nuabar  such  that  k  >  — x{^,IBa(u0)la>  ^ j ,  and  consider 
the  function 

f  »  ^^(u^)]^  _  *)  +  “  «ax{I6a(uB)]h  -  k»0)  . 

It  is  i— adiate  to  verify  that  *  9  W^' 1  (fl^)  and  therefore  it  can  be  used  as  a  test 
function  in  (7.21).  He  obtain 


/  /  ♦  v«v,h-vivv,h  - * 


(7.21) 


+  (v  u  )  ( 18  (u  ))h  -  k]+}dxdT  »  -/  [g(x,T,u  ))  (B  (u  )  -  k)+do 
a  x  a  n  nan  '  0  n  a  in 

st 


He  perform  an 
the  fact  that 


integration  by  parts  in  the  first  integral  and  let  h  ♦  0,  exploiting 
(and  hence  B^u^))  belongs  to  V*'0^)  (11).  This  gives 
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0 


(7.22) 


4  /  (8  (a  )  -  k)+2dx  ♦  //  7  K(u  )•?  (8  (u  )  -  k)+dxdT  + 

2  Q(t)  "  *  flt  X  "  *  ■  “ 

+  //  v  *7  u  (8  (u  )  -  k)+dxdT  -  -/  g(x,t,u  )(B  (u  )  -  k)+dO 
-  m  x  m  m  m  _  in  n  n 

8t  st 

♦  4/  (6  (u  (x))  -  lc> '*’2ci3t  . 

*  fl  "  ° 


By  our  choice  of  k,  the  last  integral  in  (7.22)  is  zero.  Me  treat  the  remaining  terms 

follows.  First  observe  that  by  our  construction  of  8  (•)  we  have 

m 

6  (s)  »  s  for  s  >  —  . 
a  m 

Therefore  since  k  >  —  we  have  (8  (u  )  -  k)+  »  (u  -  k)+.  From  this  and  routine 
b  n  b  m 

calculation  it  follows 

//  V  K(u  )-V  (8  (u  )  -  k)+dxdT  -  //  K*(u  )|V  (u  -  k)+|2dxdt  >  o  . 
n  x  a  x  m  a  ^  nxn 

t  t 

For  the  term  involving  the  velocity  we  have 


//  v^V^tBJuJ  -  k)+dxdT  -  //  vm-7v(um  -  k)+(Uiii  -  k)*  dxdt 


4  JJ  v  * V  [ (u  -  k)*]2dxdi  -  0 

2  fl  ■  *  ■ 

t 


since  div  v  «  0. 

m 


By  monotonicity  of  g(x,t,u|a)  at  the  origin 

g(x,t,u  )(u  -  k)+  >  0  . 

n  m 

Carrying  these  estimates  in  (7,22)  and  dropping  the  non-negative  terms  we  obtain 


/  (u  -  k)+  dx  <  0  a.e,  t  6  [0,T). 

0(t)  “ 

This  implies  u^fx.t)  <  k,  a.e.  (x,t)  8  0^.  The  bound  from  below  is  derived 


analogously.  By  taking 
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It  »  Mxll,  SUp  18 _<»«)•«}  » 

a©* 


the  constant  in  Law  7.3  is  aada  independent  of  a. 

This  rasult  will  bs  aoployad  to  prove  the  following  lsaaa 
La— a  7.4t  Let  K  ba  a  compact  of  0^.  Than  thara  exists  a  constant  c  depending  upon 
dist(K,3Q^)  and  c  but  not  upon  a,  such  that 

lv  I  «  ^  C(C|K)f  B  “  1.2.#*.  * 

*B  w  §  K 


Raaarki  la— a  7.4  says  that  ths  sequence  { v^}  is  uniformly  bounded  on  compacts  of  Q?. 

Ha  stress  the  fact  that  the  bound  does  not  depend  on  a.  but  depends  upon  the  aiae  of  the 

penalty  term  *”^H£(ub)  (and  hence  upon  e ) .  Also  there  is  no  claim  of  uniform 

boundedness  of  v  over  all  Q_,  but  only  on  compacts  K  C  Q  . 

m  T  T 

♦ 

For  each  a  €  ■  we  define  the  vorticity  of  v^  as  the  skew-sy— etric  tensor 

<m*3 )  of  entries 


u 


iJ 


a 


i.j  -  1.2  • 


Let  G,  G',  G"  be  regions  in  0  such  that  G"  C  G',  G'  CG,  G  CD,  and  consider  the 
cylindrical  domains 


Q  H  g  *  [tj.tjli  g1  S  O’  *  [t*,t2]>  g"  s  g"  *  tt^,t2J 
where  0  <  t1  <  t^  <  t“  <  tj  <  T. 

♦ 

The  following  local  representation  of  v  ,  will  play  a  role  in  what  follows . 

IB 

Lemma  7.5t  Let  x  +  C(x)  «  Cg(G)  such  that  C(x)  =  1  on  G'.  Then 


(7.23)  t(x)v  (x,t)  -  /  C(y)V  h(x  -  y)  /v  w_(y,t)dy  ♦  A  (x,t)  , 

*  a  »  »  m 

where  H(* )  is  the  fundaaentel  solution  of  the  Laplace  equation  and  for  x  e  G",  A  (•,•) 

n 

loc 

is  harmonic  in  x  and  (0,T),  uniformly  in  a. 
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% 


Proof  of  Lemma  7. Si  Assume  first  that  e  c"(0T).  Then  from  -&H(x  -  y)  ■  in 
2 

0' (R  ),  where  4^  denotes  the  Dirac  mass  concentrated  at  x,  we  obtain 

C(x)v‘l,(x,t)  -  /  V  h(x  -  y)*7 wfC(y)v‘i,<y<t))dy  - 

r  y  y  ■ 

2 

f  /  C«y)  -  y>  It-  ^’(y-tldy  ♦  /  [7h(x  -  y^Vciynv^’ty.tldy 

j-1  G  ayj  yj  G  * 

*  /  C<y>?  H(x  -  y)  (y,t)dy  ♦  A^Ux.t)  , 

r.  ”  1  ■ 


where 


A(i>(x,t)  »/  [Vh(x  -  y)*Vt(y)]v*l,(y,t)  -  /  -j — H(x  -  y)Vc(y)*v  (y,t)dy  . 

"  G 


G  #yi 


By  a  density  argument  this  representation  holds  for  €  Jjlflj).  Because  of  the  choice  of 


loc. 


the  cutoff  function  (,  A  is  harmonic  in  G*  and  1^  (0,T)  uniformly  in  m. 


Remark;  Such  a  representation  is  similar  to 


2  of  [ IS] .  The  point  here  is  to  point 


out  that  since  A  is  harmonic  in  X  €  G"  and  l£OC(0,T)  uniformly  in  m,  we  have 
m 

D®  A  e  lli00^-,)]  2,  uniformly  in  m  . 
x  tn  ■  T 

The  bounds  will  depend  on  dist(G* ,G" ) ,  the  constant  in  (7.3),  and  the  multi-index  a. 
Proof  of 


7.4:  Since  v  satisfies  (7.2),  denoting  with  v  a  mollification  of 
■  m  a,n 


■ 

(7.24) 


there  exists  a  differentiable  function 


Pmjh 


such  that 


—  v  -  vAv  . 
•  t  m,h  m,h 


-7 


x?m,h 


in  Q_  (see  (18]). 

T 

By  taking  the  curl 

(7.25)  57  -  vA«  *  curl(f(u  )  -  e  'H  (u  )v  ]  in  0  . 

o t  Bpti  ifti  ■  €  a  a  n  T 

♦  • 

He  already  know  that  u^  e  1^(8^,)  uniformly  in  m,  and  that  v^  €  uniformly 


+  l?(V  -  ‘‘VVVh 


.-1 


*  -  f<u  )  -  e 


Wvm 


in  m.  Consequently  setting 


we  have 


*  .  «  L,<0_)  uniformly  in  a  and  h, 

»,n  2  t 

and  there  exists  a  constant  C(c)  depending  upon  ess  suplu^J,  the  Llpschltz  constant 
Kj  of  f(*)  in  CXjl,  0Ti  vQ(x)  and  *  auch  that 

•  •  J,  o  *  C(e),  v  *'h 

mtn 


Conatruct  a  cutoff  function  (x,t)  *  C(x,t)  «  C  (Q) ,  auch  that  C(x,t)  SI  on  $}', 

C(x,t.)  So,  *♦  C(x,t)  e  c"(G)  and  0  <  C  <  1.  Than  u  .  C2  €  C*(fi)  and  vanishes  on 
1  0  M$n 

tha  parabolic  boundary  of  Q.  Multiplying  (7.34)  by  u  ,  ? 2  and  intagrating  over  Q, 

m,n 

aftar  atandard  calculations  we  obtain 


(7.26)  ess  sup  la  Cl*  {  }  ♦ 

t1<t<t2  *,h 


I  |V  ai 

x  a 


,h*5*  2,0 


C(,Vh*2.Q„ 


'"■.h1 2,Bt> 


for  a  constant  C,  dapanding  upon  c ,  diat (Q,Q' ) ,  tha  data  but  indapandant  of  a  and 
♦  ** 

h.  8inca  vb  «  J ^ (QT>  uniformly  in  a  we  have  that  a^  h  €  L2(flT>  uniformly  in  a 
and  h.  Tharsfora  from  (7.36)  racalling  tha  dafinition  of  tha  cutoff  function  C ,  we 
daduca  that  thara  axiata  a  conatant  C(C)  dapanding  upon  t  and  tha  data  but  independent 
of  a  and  h  auch  that 


(7.27) 


la 


•,h'v2(Q) 


<  C(c)  . 


Corollary  (3.2),  with  p  ■  2  iapliea  that  a  e  L.(QI>)  uniformly  in  a  and  h. 

■ffl  4 

Therefore  ainca  tha  choicaa  of  Q,  Q',  Q*  ara  arbitrary,  wa  deduce  that  e  l2oc(  $^) , 
uniformly  in  a,  with  bounds  dapanding  upon  e. 

From  tha  repraaantation  (7.23)  and  the  Calderon-Zygmund  theory  of  singular  integrals 
[20],  wa  deduce  that 

^  \nc 

V  v  e  L,  (!)_)  uniformly  in  o  , 
x  a  4  T 


-29' 


♦  1  Qq 

and  therefore  v  e  V,  (0  )  uniformly  in  m.  By  Oorollary  3.2,  with  p  •  4  we  find 

W  T 

♦  loc 

€  Lg  (Q^)  uniformly  in  m.  Mow  letting  h  ♦  0  in  P'(Q^)  in  (7.25)  we  have 

(7.28)  f-w  -viu  -  curll*  )  in  P*(n_)  . 

ot  ra  a  a  T 

loc 

Prom  (7.27)  we  see  that  m  is  a  solution  of  (7.28)  which  belongs  to  V.  (Q_).  Since 

a  2  T 

Xoc 

6  Lg  (flT>  uniformly  in  m,  standard  parabolic  theory  [11,  16,  17),  implies  that 

Xoc 

<  ly  (0^)  uniformly  in  m,  with  local  bounds  depending  upon  e. 

The  lemma  is  now  a  consequence  of  the  representation  (7.23). 

3  Xoc 

Finally  we  employ  Leans  7.4  to  show  that  e  L^  (0^) . 

Lemma  7.6i  Let  K  be  a  compact  of  0^.  Then  there  exists  a  constant  C  depending  upon 
dist(  K,3Qt) ,  e  and  the  data,  but  independent  of  m  such  that 


'It  tt.*L2((0  *  C  * 

Proof  of  Lemma  7.6i  Let  tjC  9  be  compacts  and  consider  the  cylindrical  domains 

Q1  S  x  [t^tjjf  02  2  K2  *  ftJ'ty  wh«re 

0  <  t'  <  t^  <  tj  <  tj  <  T  . 

88 

Construct  a  cutoff  function  v>( *,t )  €  C  (Q^)  such  that 
( i )  aupp  *  COj 
(ii)  *>  =  1,  ( x, t )  €  Q1  . 

Next  consider  identity  (7.21)  and  set 


K(ua>  "  *» 


and 

YU)  -  6  (*”'(»))  . 

ram  m  n» 

It  is  clear  that  by  virtue  of  the  assumptions  on  K(*)  it  will  be  sufficient  to  prove 


'It  *»'  2,  K 


<  c  . 
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In  (7.21)  with  the  indicated  change  of  variable,  choose  the  teat  function 


Since 


(7.29) 


supp  1 1  C  Qj,  the  term  involving  integrations  on  ST  drops  out,  and  we  have 

fj  tfe  It  1%V2  +  Wh  It  w*/ + 

B2 

+  7  I*  k  4r  t*  lK*v  *2  *  [v  -v  k-1(*  )).  4r  t*  JK*2}dxdt  -  o 
x  ah  3t  a  h  x  ax  a  h  3t  ah  1 


Consider  the  first  Integrand  in  (7.29)  and  recall  that  by  virtue  of  assumption  [hj] 
and  the  construction  of  the  sequence  B  ( • )  we  have 

B 

Y^(s)  -  B^Ik'1  (a)]*!!”1 '  (s)  >  X"1  s  S  r\{0}  . 

Mow  from  the  definition  of  Steklov  averaging  it  follows  that 


a  a  ♦  h>>  *  r.(t(t))  a  Jt  +  h)  -  z  (t) 

k  'W'hfc  t'.'h-  -  - - ‘ 


z  (t  +  h)  -  s  (t)  ,  . 

- n - =— 1  -S’llt^hl- 


For  the  term  involving  gradients  we  have 


U  Wh  k  Wh*2^ 

*2 


iff  k  iWh 

W«j 


2<P2dxdt 
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These  remarks  in  (7.29)  give 


!f  lit  ‘Vh1  Vdxdt  -i  //  IWJ2  It"  dXdt  ■  2  U  Wh  It  l'«W 


(7.30) 


-  U  ‘VV^V^It  tz»1h^2dxdt  “  X1  +  X2  +  *3 

«*» 


Ha  estimate  the  1^,  1  -  1,2,3  separately.  Proa  (7.3)  it  follows  that  the 

integral  I1  is  uniformly  bounded  with  respect  to  ■  and  that  the  bound  will  depend  on 
i.e.  the  distance  between  Q,  and  Q2. 

2  12 

As  for  l2  we  use  the  Cauchy  inequality  ab  <  na  +  —  b  n  >  0,  to  obtain 


|i2l  «  »  //  l|t  UB]h|V  +$fl  l%l*.1hlal%*l*«t 


2|»  |2j 

X * 


He  do  the  saae  with  I. 


IJ3I  <  1  //  lf^  [*|(]h|2«2dxdt  +^IJ  ll^M,7xK”1(,B>lh|2«2 


2  “2 
Carrying  these  estimates  in  (7.30)  we  obtain 


l^;1  -  2n]  //  (.Jh|  Vdxdt  <  c(n)  II  |Vxl«Jh|2Uktl  ♦  lVl2]dxdt  + 

®2  °2 

i«%v“,«-->ihia  2  • 


since  we  have  previously  shown  that  v_  is  locally  bounded  independent  of  a  the  last 
integral  is  uniforaly  bounded  with  respect  to  a.  Finally  we  choose  n  -  X~'/4  and  recall 
that  <t  =  1  on  Q,  •  This  yields 

'It  ‘Vh'a.Q,  <  C(,V.,Q2'di't<«l'«2)',VxV2.aT>  * 

Since  this  independent  of  h,  from  a  result  of  (11 J  it  follows  that  the  weak 

derivative  exist  and  is  a  locally  square  suaaable  function  in  nT.  The  lemma  is 


proved 


I 


Remark:  We  stress  the  following  two  facts 
3u 

(a)  fl  Lj  only  locally,  uniformly  in  m. 

(b]  The  a-priori  bound 


3u 


»*t*2.Qt  <  C(,V-,QJ'di>t  (K'3QT)',V«,2.Qt) 


does  not  depend  upon  ra,  since  '  Iv  I  and  IV  u  I,  _  are  bounded  independently 

i®  '*2  x  * 

of  a.  but  does  depend  upon  e  via  I*  L  .  . 

»  •'*2 

7-(b)  The  limit  as  ■•»•** 

We  now  conclude  the  proof  of  Proposition  6.1  by  letting  m  ♦  •  in  (7.1  )-(7.2).  Prom 
(7.3)  it  follows  that  the  sequences  (u^li  {  v^}  are  weakly  compact  in  and 

j)(Qt>  respectively!  hence  subsequences  can  be  selected  and  relabeled  with  m,  such 
that 


u  ♦  u_  weakly  in  w',0(Q_) 
a  e  2  T 

v^  ♦  ve  weakly  in  J1 (0T>  . 


Lemma  7.1 1»  Let  e  >  0  be  fixed.  There  exists  a  subsequence  (relabeled  with  m)  such 
that 

u^  ♦  u£  strongly  in  L^O^)  . 

Proof  of  Lemma  7.11. 

Let  K  be  a  compact  of  nT.  Then  by  Lemma  (7.3)  and  (7.11)  we  have 

K  +  '7TI2,(C+,VxV2.IC<c<K> 

where  C  depends  upon  e,  dlat(K,dS)^)  but  not  upon  m. 

Therefore  for  a  subsequence 

u  ♦  u„  strongly  in  L,(K)r  and 
m  c  c 

um  *  “e  a,a*  tn  K  • 

Now  by  the  uniqueness  of  the  weak  limit  ue  we  have 

ue  =  u£f  a.e.  (x,t)  8  K  . 


* 
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w 


Let  K  ■  sequence  of  compacts  of  Q„  such  that  K  C  K  and  u  K 
*  T  p  p+-1  1 


p>1 


Qt  .  Than  by  a 


diagonalication  process  a  subsequence  can  be  selected  and  relabeled  with  a  such  that 


u£  a.e.  In 


Since  lu  l—  -  is  equibounded,  by  the  Le basque  doeinated  convergence  theoree  we  have 
'  T 


u  strongly  in  LjfSl^)  * 


By  the  trace  theoree  and  (7.3) 

trace  u^  *  trace  u£  strongly  in  L^tS^)  ' 

and  by  monotonicity  of  8 < • » ,  S^tu^)  ♦  w  8  B(u)  weakly  in  L2(QT). 

be— a  7.12t  For  a  subsequence  (relabeled  with  n) 

♦  ♦ 

v^  ♦  v£  strongly  in  J(QT>  . 

Proof  of  La— a  7.12.  A  consequence  of  Le— a  7.  It  is  that  f(u  )  ♦  f(u  )  and 
"  B  E 

■gtu^)  ♦  H£ ( u£ )  strongly  in  b2(flT).  Therefore  since  the  space  dimension  is  N  -  2  and 

*  >  0  is  fixed,  the  strong  convergence  of  va  to  v£  follows  by  a  straightforward 

adaptation  of  Serrin’s  stability  theorem  (Theorem  6  of  [19]  page  S3). 

We  can  now  pass  to  the  limit  in  (7.1),  (7.2)  as  a  ♦  »  (for  m  labeling  the 

particular  subsequence  chosen  above)  to  obtain  the  existence  of  a  pair  (u  ,v  )  such  that 

c  c 


(7.31) 


V  (Q  )' 
2l  T' 


■V  - 


<  c„ 


vv 


where  the  constant  CQ  does  not  depend  on  e,  and  u£,  v£  satisfy 

(7.32)  //  {-B(ue)v»t  +  VXK(UC>‘VX*  +  Vg*VxUg*}dxdt  -  -/  g(x,t,Ug)*dc  +  /  8(u0)*(x,0)dx 

°T  ST  0 

for  all  *  8  w2'1(HT)  such  that  *(x,T)  -  0,  and 

//  +  W  v  i  V  I  +  e-1H  (u  )v  »$)dxdt  - 

Q  »  t  »  W  X  t  t  b 

T 

“  //  f(u  )$  dxdt  +  /  v  (x)$(x,0)dx 

6  fl  0 

T 

for  all  |  hx.D  s  0. 
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Statements  (b)-(e)  are  obvious  from  the  estimates  we  have  established  in  Lemmas  7.3, 
7.4,  7.10,  7.2. 


8.  PROOF  OF  PROPOSITION  6.2 

Since  the  modulus  of  continuity  has  to  be  uniform  in  e 
estimates  that  we  have  which  are  uniform  with  respect  to  e. 


start  by  listing  the 


'V-  fl  <  c2'  V  e  >  0  . 

'  T 


,ue'v2<nT)'  'Vv2<aT>  *  co  * 


From  the  second  of  (8.2)  and  Corollary  3.1  we  have 


lv  I  -  <  (Const  Independent  of  e)  . 

E  4*"t 


Also  the  qualitative  information  in  (d)  in  Proposition  6.1  is  essential  in  order  to  justify 
soma  of  the  calculations  below,  but  the  modulus  of  continuity  will  not  depend  on  the  local 
estimates  for  '§*  uf'2,K* 

The  function  u  (•)  claimed  by  Proposition  6.2  for  will  only  depend  on  the 

quantities  listed  in  (8. 1 )— (8.3). 

In  this  section  the  restriction  N  ■  2  will  play  a  dramatic  role.  The  flow  of  the 
proof  is  like  the  arguments  produced  in  (7,  8].  Now  the  order  of  summability  (8.3)  of  v£ 
is  not  high  enough  to  fulfill  the  assumptions  of  Theorem  1  of  (7).  For  this  reason  a 
modification  in  the  proof  is  needed  where  we  will  exploit  both  the  dimension  N  ■  2  and 
the  particular  structure  of  the  equation  corresponding  to  identity  (6.1). 

Since  the  arguments  have  been  presented  in  detail  in  (7)  we  will  limit  ourselves  to 
pointing  out  the  differences  that  occur  at  various  steps  in  the  proof. 

The  main  idea  will  consist  in  showing  that  given  (xg*tg)  8  we  can  construct  a 

sequence  of  cylinders  ^  •centered"  at  <xQ, tQ),  such  that  0  <^+1  and  shrinking  to 
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U0,t0),  where  the  oscillation  of  u£  decreases,  according  to  the  operator  associated 
with  (6.1)  and  in  a  way  deterained  only  by  the  quantities  in  (8.1)-(8.3).  Given  n  >  0, 
this  process  will  also  prescribe  the  size  of  a  cylinder  Q(n),  where 


ess  osc  u  <  n,  V  e  >  0  . 

C(n)  e 


This  will  obviously  yield  a  modulus  of  continuity  for  u£  over  coapacts  K  C  0T,  in  an 
uni fora  fashion  with  respect  to  £. 

[8. A]  Preliminary  material. 

First  we  report  a  result  of  [11]  which  in  our  context  can  be  stated  as  follows.  Let 
u  satisfy  (6.1)  and  u  8  (fl  ).  Then  u  satisfies  the  integral  identity 

E  £  a | IOC  T  £ 


2  2  + 

(8.4)  /  6(ue)^(x,t)dx  ♦  /  /  {-S(u£)^t  ♦  7xK(ue)‘7x^  ♦  Vevx“£^}dxdt 


ti  a 


for  all  <fi  9  W^'1^)  auch  that  x  ♦  *(x,t)  is  supported  in  Q  for  all  t  8  [t^.tjl, 
and  all  intervals  (t^tjl  C  (0,T).  We  will  consider  cylinders  contained  in  of  the 

following  special  fora.  Let  ( xQ , tQ )  be  an  arbitrary  point  in  0^  and  denote  with 
B(R)  the  ball  {|x  -  xQ|  <  r)  and  with  £^(8)  the  cylinder 

5^(8)  =  B(R)  x  (t()  -  8R2,tQ]  . 

Also  if  °1»°2  8  (0,1)  we  set 

Ve.°r° 2>  5  B(R  "  °iR)  *  lto  "  8(1  "  <,2>R2,to)  • 

Consider  the  definition  (4.2)  of  the  graph  f)(*)  and  set 
(8.5)  0(u£)  -  u£  +  H(u£) 

where  s  ♦  H(s)  is  the  graph 


(8.6) 


H(s) 


0 

(-L,0]  , 
-L 


S  >  0 

s  a  0 

S  <  0 


In  (3.4)  we  employ  a  test  function  <i(x,t)  supported  in  the  ball  B(R)  for  all 
2 

t  8  (tg  -  8r  ,tg]  where  R  and  8  are  assumed  to  be  so  small  that  0^(8)  C  By  the 

3 

results  of  the  previous  section  u£  9  L2(QR(8)),  therefore  substituting  (8.5)  in  (8.4) 
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and  integrating  by  parts  with  the  indicated  choice  of  V  we  obtain  the  identity 


(8.7) 


,t  t 

/  H(u  )<P(x,t)  |  ,  -  /  ,  /  H(u  )*  dxdt  + 

8(R)  tQ-0R  t0-8R  B(R> 


+  /  ,  /  ilrV  +  V vK(ue,*Vx  +  V  *v  u  ^Jdxdt  =  0 

t0-eR2  b(r)  3t  e  *  6  x  e  x  e 

for  all  *>  6  * 1  £  Qr  f  8 ) )  *  and  all  t  8  ItQ  -  8R2,t0).  The  purpose  of  (8.7)  is  to  isolate 

the  contribution  coming  from  the  jump  in  8(»)  with  respect  to  the  rest  of  the  equation. 
Next  we  construct  particular  te3t  functions  in  (8.7). 

Let  (x,t)  *  C(x,t)  be  a  cutoff  function  in  QR(0)  satisfying 


(i)  C(«,t)  ec0(B(R)),  |vxcj  <  (olR) 


-1 


(ii)  C(x,tQ  -  0R2)  =0,  x  e  B(R) »  0  <  <  (0o2R2)-1, 

(iii)  C(x,t)  5  1,  (x,t)  e  Q^.o,,^). 

Let  k  6  R  and  consider  the  functions 


(ue  -  k)  +  »  max{ue  -  kjO) 


(ue  -  k)  »  raax{-(u£  -  k);0} 


It  is  obvious  that  if  u  G  L  (t(0)),  then  (u  -  k)*  e  L  ( Q  ( © ) ) ,  r,s  >  1.  It  is 
c  r  f  s  R  t  1 1 3  R 

known  that  if  u  6  1 (Q  (8 ) )  then  also  (u  -  k)*  6  w' ' 1 (Q  (8 ) ) ,  (cf.  llll). 

w  4  K  t  4  R 


In  (8.7)  we  will  choose 


*(x,t)  -  ±(u£  -  k)*CZ(x,t) 


For  simplicity  of  notation  we  will  drop  the  subscript  e  and  set 


-**(k,t  >  -  /  ±H(u) (u  -  k)±t2(x,t) 


B(R) 


V0R 


/  ,/  H(u)(±(u  -  k)±C2(x,r)]  dxdT  . 

tQ-0R  B(R)  * 
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The  regaining  teres  in  (8.7)  are  transformed  as  follows 


f  j  ±  *  uOl-kjVu.Ttfxdt  -l/t  /  J-  [<U-k)t)V(x.T)dXdT  - 

tQ-8R2  BOO  3t  2  t0-9R  B(R) 


-4-/  I(u-k)±)2c2(x,t)dx  -  4  /  j  2  /  [(u-k)*]2C(x,T)  C  dxdt  > 

2  B(R)  2  tQ-8R  B(R) 


>l|(u-)0±CI2<B(R)(t,  -  I  (u-k)*(Ct)^  *2,^(8)  * 


For  the  tere  involving  V  ^  we  have 


„  J  ±V  K(u)V  (u-k)*C2dxdT  -  /  /  K*(u)|V  (u-k)*!2?2  + 

to-0R2  BCR)  X  X  to-0R^  B(R) 


+  f  ,  /  2K’(u)V  (u-k)*C(u-k)*V  C  >  X  /  J  |V  (u-k)t|2C2dxdT  - 

to-0R2  B(R)  X  *  to-0R2  B(R) 


-e/  J  |V  (u-k)±|2C2dxdT 

tQ-0R^  B(R) 


//  Uu-k)*]2|V  C|2dxdt  . 

Qr(0) 


(*0  “  2 
to-0R2 


/  |V  (u-k)*|VdxdT  - 

B(R) 


l(u-k)t|VxC|l 


2,Qr(6> 


Me  treat  the  term  involving  the  velocity  as  follows 


r 

\ 


I  .  /  ±v*V  al {u-k)*K2(x,T)dxdT  » 
tQ-6R  B(R)  * 


if  ,  J  £•(?  l(u-k)*)2)C2(x,T)dxdT 

•  -  An*  n/h»  X 


t0-6R  B(R) 


-  4  /  ,  2  /  $[<U-*)4J2CV  C  dxdt  . 

tQ-8R  B(R) 


Confining  these  estimates  as  parts  of  (B.7)  we  obtain 


\  I (u-k)*CI2  ,<t)  ♦  (X  -  0  /  2  /  |V  (u-k)4|2C2dxdt  < 

2  2,B'R)  0  t0-8R2  b(r)  x 


•ve> 


ve> 


+  /  ,/  |v||V  C||<U-*)*|2dxdT  ♦  *t(k,tn) 


to-0R  B(R) 


for  all  t  «  [tQ  -  8r  ,tQJ. 


Now 


choose  e  “  "2  recall  the  construction  of  the  cutoff  function  ? (x,t)  and  use 


the  arbitrarity  of  t  8  ItQ  -  0R, tQl ,  to  conclude  that  there  exists  a  constant  y 


depending  only  upon  the  data  such  that 

l(u-k)*l2  „  <  Yl(o.R)*2  +  (0  9r2 )”^ 1 • I (u-k )41 2  + 

V2  1  2  2,fiR<9) 

(8.8) 


♦  Y(o,R)*’  //  jv(((u-k)*?2dxdr  +  sup  y^Ck.t.) 

Qr(8,  «"V8r2,V 
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Inequalities  (8.8)  are  valid  Cor  every  cylinder  QR(8)  c  ^T>  every  pair  a^.c^  6  (0.1) 
and  every  real  number  k.  They  are  one  of  the  main  tools  in  the  proof  of  continuity. 
Another  tool  is  the  following  logarithmic  estimate. 

(8.B]  A  logarithmic  estimate. 

Lemma  8. 1 :  Let  k  8  R*.  u  >  ess  sup  (u  -  k)+  and  0  <  n  <  (i.  Set 

V9) 

♦(x,t)  «  tn+  J- - - - — - 1  «  max  /  In  - - - — - 1  oj  . 

|m  -  (u  -  k)  +  nj  [  u  -  (u  -  k)+  +  n  J 

2 

Then  there  exists  a  constant  C  -  C(8)  such  that  for  all  t  8  [tQ  -  0R  ,tQ) 

/  *2(x,t)dx  <  /  ^(x.t  -  8r2)  +  ^-  (tn  £)  meas  B(R) 

B(R-V>  B(R)  *? 

Remark »  For  simplicity  of  notation  we  will  use  the  same  symbol  8  for  t(x.t)  and 

•<w  i 

♦(u(x,t)).  in  what  follows  ♦  '  will  mean  ♦  . 

In  the  cylinder  QR(8 )  construct  a  cutoff  function  x  ♦  C(x)  such  that 
(i)  C(x)  Cq(B(r)  ) ,  |Vxc|  <  (o^)-1 
(ii)  C(x)  SI,  x  8  B(R  -  0tR). 

Proof  of  Lemma  8.1.  In  (8.7)  consider  the  following  test  function 

<P  (xf  t)  »  (t2)'?2!*) 

where  x  ♦  C(x)  is  as  above. 

It  is  apparent  that  *p  6  W^'tO^B))  and  that  (f2)*  ■  2(1  +  <>)W)2. 

2  + 

Since  (♦  )'  vanishes  at  those  points  (x,t)  8  (^(8)  where  (u  -  k)  <  n  and 

n  >  0,  the  terms  in  (8.7)  involving  H(ue)  do  not  give  any  contribution.  The  term 
du 

Involving  gives 

/  -/  u(t2)'c2(x)dxdt  -  /  t2(x,T)c2(x)dx  I  . 

t0-®R  B(R)  B(R)  to"®* 

We  estimate  the  remaining  terms  as  follows 
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/*  /  K'(u)V  ut (2d  ♦  ♦){*,)2}?uC2  ♦  <*2),V2jdxdt  * 

tg-flR2  8(10  X 


>  1  f  ,  /  2(1  ♦  +)<♦• )2|V  u|V  +  J 

0  t0-6R2  B(R) 


where 


J 


4  _  /  KMu>*  ♦•V  u  CV  C  dxdT  <  2*  /  ,/  ( 1  +  7) (*’ )2| 7^1 2? 

t0-«R2  B(R)  X  X  t0-6R2  B(R) 


■f 


/ 

v9R 


;  »|7  c|2dxdt 

B(R) 


For  the  tern  involving  the  velocities 


we  have 


2  J1  „  /  S*7  ut'tC^XdT  -  2  /  ,  /  v»V  ♦•♦C2dxdT  < 

t0-BR2  B(R)  x  t0-«R2  B(R) 


«  2c  /  /  (1  ♦  ♦)(♦') 


t0-8R‘ 


B(R) 


l|Vxu|2C2dxdT 


j  7|  v|  2C2dxdt 
t0-8RX  B(R) 


i( 


Collecting  these  estleates  we  have 


•*d 


2 

2,B(R) 


(t)  ♦  (2  X  -  4*)  /  , 

V8R 


/  (1  +  *)(*')2|7  u|  VdxdT  < 

B(R) 


<  l«l 


2 

2,B(R) 


(tQ  -  0r2> 


air  , 

*T2  t.-SR2 
0  0 


/  7(7  C|2dxdt  ♦ 

B(R) 


/  ,/  t|v|Vdxdx 

q“8r‘  B(R) 


(8.9) 


Now  choose  e  «  and  observe  that  ♦  <  In  Recalling  the  construction  at  C  we 

4  n 

conclude  that  there  exists  a  constant  C  depending  only  upon  the  data  suah  that 

/  *2(x,t)dx  <  /  *2(x,t  -  8R2)dx  ♦ 

B(R-OjR)  B(R) 

(8.10) 

+  C  in  jj-  { (<>1R)”28R2  meas  B(R)  +  lvel2  g  [8r2  aeas  B(R)  1^*  }  . 

+  2 
Since  lv_l  .  is  unifornly  bounded  with  respect  to  c  and  [8r  aeas  B(R))/*  < 

e  4«“t 

<  (const)  aeas  B(R),  for  N  »  2,  it  follows  fro*  (8. 10)  that. 


/  *2(x,t)dx  <  /  *2(x,t  -  8R2)dx  ♦%(*»£) 

B(R-OjR)  B(R) 


Mas  B(R)  . 


The  lenna  is  proved. 

Let  us  return  now  to  the  inequalities  (8.8)  and  estinate  the  tan  involving  velocities 
as  follows. 


„(t)  =  (x  6  B(R)|(u(x,t)  -  k)*  >  0} 


M(k,R)  «  ess  sup  (u  -  k)1 


Tf(o.R)  If  |v|((u  -  k)VdxdT  < 
0,(8) 


<  Y(a.R)'1lM(k,R)l2  //  | v|xt (u  -  k)*  >  OJdxdt 

v9> 

where  x((u  -  k)*  >  0]  is  the  characteristic  function  of  the  set  l(u  -  k)*  >  0)  n  {^(8). 


tQ 

//  |v|xl(u  -  k)*  >  0)dxdT  <  ijl  a  •  /  ,  nsas  _(t)  ,/4 

qr(8)  4'8t  [t0-8R2  *'R 
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Therefore,  by  changing  the  constant  7  appropriately,  (8.8)  can  be  rewritten  as 


Ku  -  JO*!2 


,  q  <  r[«y.R)~2  +  <a,9R2)"Vl<u  -  k)*l2  ...  + 

V2  lQR(9'Cl'a2))  R 


<8.12) 


+  Y^R)"’  lK(k,R)]2 


/  meas  A*  _<x)dx 

LtQ-8R2  k,R 


3/4 


sup  T*  (It,  t_>  . 

t«tv6R  'V 


Let  us  now  show  how  to  conclude  the  proof  of  Proposition  6.2. 

In  [7.8]  we  demonstrated  that  the  existence  of  a  modulus  of  continuity  for  weak 
solutions  of  singular  quasilinear  parabolic  equations  in  divergence  form  was  solely  a 
consequence  of  inequalities  (2.7)  page  16  of  17]  and  lemma  2.2  page  19  of  [7].  Mow  the 
analog  of  Lemma  2.2  of  [7]  is  precisely  Lemma  8.1  here.  Actually  the  structure  of  our 
equation  leads  to  a  less  complicated  logarithmic  estimate.  We  stress  the  fact  that  the 
derivation  of  Lemma  8.1  employs  in  an  essential  way  the  fact  that  the  dimension  N  is  2. 

As  for  inequalities  (2.7)  of  (71,  their  analog  here  are  inequalities  (8.12).  There  is 
only  a  slight  difference  in  the  term  involving  A*  R(t). 

In  [71  such  term  reads 


I  - 


r  ±  q 

/  ,  [meas  A*  D(T)pdT 


V8R 


k,R 


7  (1+0 


where  q,r  >  0  are  linked  by 

IN  N 
—  +  —  »  — 
r  2q  4 

and  r  8  (0,1 ). 

Therefore  I  can  be  estimated  by 
(8.13)  I  «  Const  RN*RNlt 

In  our  case  the  analogous  term 


i*  -  toT1*”1 


/  meas  AT  „(t)dt 


V 


3/4 


iM(k,R)r 
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is  estiaated  aa 

(8.14)  I'  <  const  (M(k,R)]  V  .  (8*2)  • 

Mot a  that  in  this  astiaate  too  N  *  2  is  essential. 

Finally  let  us  show  the  difference  between  (8.13)  and  (8.14)  does  not  affect  the  proof 
in  (7J. 

Let 

*  ess  sup  u;  y  *  ess  inf  u 

aR(0)  V6) 


so  that  the  oscillation  u  of  u  in  (^(9)  is 

u  *  u*  -  u"  . 

Inequalities  (8.12)  are  eaployed  with  the  choice  of  k  given  by 

.  +u) 


.  -  a) 

k  *  u  ♦  — 


where  s  0  8.  Consequently,  froa  the  definition  of  (u  -  k)  , 


M(k,R)  <  —  . 
2 


(l)  /  2 

Now  if  —  <  R  '  ,  the  oscillation  can  be  bounded  with  a  power  of  R,  and  there  is 
28 

nothing  to  prove.  The  case  to  exaaine  is  then  when 

—  >  rn*/2 
28 

r  u)  \  2 

In  [7]  we  estimated  R  from  above  with  [ — J  and  carried  out  the  arguments  with  such 

2® 

an  estimate.  Now  this  is  precisely  the  content  of  (8.14)  via  (8.15).  The  tera 

(M(k,R))2  in  (8.14)  therefore  plays  the  role  of  R***  in  (7)  when  division  by  (— )2  is 

2* 

carried  out. 

We  omit  the  details  (which  are  given  in  l 7 ] )  noticing  that  our  situation  is  in  fact 
easier  due  to  the  simpler  structure  of  the  equation.  The  proof  is  complete. 
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ggMrti  The  assumptions  uQ  €  L^SUO))  and  *C  ♦  g(x,t,£)  monotone  at  the  origin*,  were 
used  in  the  proof  of  lemma  7,3,  Now  it  is  apparent  from  the  previous  arguments  that  one 
only  needs 

Ug  e  i^oc(aT)  independent  of  m  and  c  . 

The  latter  can  be  proved  starting  from  inequalities  (8,8)  with  the  aid  of  Theorem  6.2  of 
[11]  page  105.  Consequently  one  needs  only  to  assume  uQ  €  L2<Q(0))  and  the  monotonicity 
condition  on  g(x,t,* )  at  zero  can  be  relaxed. 
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